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SUMMARY 


The occurence of unstable vibrations in rotating machinery requires 
corrective “^asures for improvements of the stability behavior. Before a measure 
will be realized, different possibilities should be investigated. In this paper 
a simple approximate method is represented, to find out the influence of para- 
meter changes to the stability behavior. The method is based on an expansion of 
e eigenvalues in terms of systemparameters . Influence coefficients show in a 
very intuitive way the effect of structural modifications. The method first of 
ail was applicated to simple nonconservative rotor models. Furthermore it was 
approved for an unsymmetric rotor of a test rig. 


INTRODUCTION 


. occurence of instability in rotating machinery may be caused by 

different effects, for example oil film forces in journal bearings, forces in 
sea s, unsymmetric shafts, internal damping etc. A machine designer wants to 
know whether a rotor will run stable during operation and what size the stabi- 
ity threshold speed will have. Furthermore he needs information about the pa- 
rameters influencing the stability behavior. 


Important informations about the stability of a linear rotor system can be 
obtained by calculations in the design stage. Because of uncertain input data 
or calculation, the results have to be considered critically. It is possible 
that unstable vibrations may occur during operation, although the calculation’ 
was predicting a stable machine. On the other hand parameters may change during 
operation leading to increasing oscillations. In such cases suitable corrective 
measures for improvements of the stability behavior are required. Before a 
measure will be realized, different possibilities should be investigated finding 
out the simplest and most effective one. For that it would be very useful to 
have approximate formulas, expressing the sensitivity of the dynamic behavior to 
changes of system parameters. 


In linear rotor systems with nonconservative effects the stability can be 
valuated by means of the system eigenvalues. The real parts of the eigenvalues 
etermine, whether the natural motion is decreasing or increasing. If the 
variations of eigenvalues caused by variations of systemparameters (mass-, 
amping , stiffness coefficients) are known, an estimation of system modifi- 
cations to the stability behavior is possible. Such sensitivities, respectively 
influence coefficients, expressing the change of eigenvalues to changes of 
systemparameters are presented here. 
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Lund (ref. 1) has developed a method to calculate sensitivities of the 
critical speeds (eigenfrequencies) of a conservative rotor to changes in the 
design. Dresig (ref. 2) gives a more general development for conservative mecha- 
nical structures. The basic idea in his method is an expansion of the eigen- 
values in terms of the system parameters. In Taylor's expansion derivatives of 
the eigenvalues to the system parameters are needed. Such derivatives were 
developed from Plaut and Huseyin (ref. 3) and from Fox and Kapoor (ref. A). 

Based on Taylor's expansion for complex eigenvalues in this paper an 
improvement of the method is presented for rotor systems with nonconservative 
mechanisms (oil film bearings, seals ets.) Linear, quadratic or higher order 
formulas are obtained, depending on the order of the derivatives taken into 
consideration in the expansion. With the linear formula very simple influence 
coefficients can be defined, pointing out the influence of special parameters 
to an eigenvalue. Superposition of different parameter changes is possible in 
this special case. Using the formulas the eigenvalues, the left-hand and 
right-hand eigenvectors of the original system (without parameter changes) 
must be known. 

The application of the method is demonstrated for simple nonconservative 
rotor models, investigating the effects of stiffness and damping coefficients 
to the stability behavior. Furthermore the influence of the mass variation to 
the eigenvalues of an unsymmetric rotor is determined. The predicted vibration 
behavior caused by parameter variations could be confirmed by measurements. 


NATURAL VIBRATIONS OF LINEAR ROTORS 


A turbomachine consisting of a high pressure turbine, a generator and an 
exciter is shown in figure 1. The power of the machine is 110 MW and the opera- 
ting speed 3000 rpm. Contrary to nonrotating structures the dynamic behavior of 
such rotating machines is influenced by additional effects. Of great importance 
are selfexciting and damping effects (nonconservative effects), caused by the 
oil film forces of journal bearings, forces in seals etc. 

In linear rotor dynamics the natural vibrations can be described by linear 
equations of motion, usually derived by means of the finite element method. 


Mq + Cq + Kq = 0 ( 1 ) 

M mass matrix (order NxN) 

C damping matrix (order NxN) 

K stiffness matrix (order NxN) 

q displacement vector 
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The equations express the equilibrium of inertia, damping, and stiffness forces. 
Because of the nonconservative effects damping and stiffness matrix contain also 
skewsymmetric and nonsymmetric terms besides the symmetric ones. Furthermore some 

of the matrix elements depend on the operating conditions of the machine (speed 
power, pressure etc.). ^ ’ 


Investigating the natural motion (stability) of 
to be solved. The solution has the form 


a rotor equation (I) has 


q(t) = q e^*" 


( 2 ) 


Substitution yields the quadratic eigenvalue problem 
2 

{A M + Ac + k) q = 0 


(3) 


with 2N eigenvalues A^ and corresponding eigenvectors (natural modes) q . The 

eigenvalues as well as the eigenvectors mainly occur in conjugate complex 
pairs 


eigenvalues: 

A = a + ioj 
n n n 

, A = a - 

n n 

io) 

n 

(4) 

eigenvectors : 

q = s + it 
-n -n -n 

-n — n 

it 

-n 

(5) 


The part of the solution, which belongs to such a conjugate complex pair, can 
be written as r t- > 


^ at 


( 6 ) 


w IS the circular natural frequency and a the damping constant (decay constant) 
The damping constant, ^respectively the rea? part of the eigenvalue determines, 
whether the solution q^(t) decreases (a^ < 0 ) or increases (a > 0 ). 


f lowest natural frequencies and the corresponding damping constants 

o the turbomachine (figure 1 ) are plotted versus the running speed in figure 2 . 
e eigenvalues are changing with speed, especially the damping constants a, and 
3 . he diagram shows that the rotor instability onset speed is 3400 rpm deter- 
mined by the zero value of aj. The relative distance between the instability 
onset speed and the operating speed is very low. Therefore unstable vibrations 
may occur if additional destabilizing forces are acting on the rotor. If such 
calculated results are known in a machines design stage, corrective measures 
improving the stability behavior should be arranged. If instability occurs in 
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operating machines similar problems have to be solved. Approximate formulas, 
respectively influence coefficients expressing the influence of parameter 
changes to the eigenvalues (stability behavior) may be very useful for the 
above mentioned requirements. 


INFLUENCE OF PARAMETER CHANGES TO THE 
STABILITY BEHAVIOR 


Eigenvalue Problem and Modal Parameters 

Our aim is to find formulas for an approximate calculation of eigenvalue 
variations caused by system modifications (parameter changes). Besides the 
designed parameter changes this formulas will contain modal parameters of the 
initial dynamic system, that means eigenvalues A, left-hand eigenvectors ^ and 
right-hand eigenvectors q of the system without parameter changes. Normally they 
are known from calculations (ref. 5). 

Complementary to the eigenvalue problem (eq. (3)) the corresponding left- 
hand-eigenvalue problem with transposed matrices 


{xV + + K^} 1 = 0 


(7) 


has the same eigenvalues A, but the mentioned left-hand eigenvectors il. Working 
with the procedure of inverse vector iteration for complex eigenvalues (ref. 5) 
a desired number of the modal parameters A, q, can be calculated effectively. 
For further derivations we suppose that the interesting modal parameters of the 
initial system are given. 


Taylor’s Expansion for the Eigenvalues 

The basic idea of the approximate method is an expansion (Taylor's series) 
of the eigenvalues in terms of the generalized system parameters Pj^ 

(k = 1,2, ...K) 


X = X 
n no 


+ E 


i{^ 

r " 


Ap j + 


9p. 


^^2 3p, 


Apj^}’^ ^(P] . V>2*‘ 


( 8 ) 


in which the p, may be mass, damping, or stiffness parameters.^ 

X is the changed eigenvalue after a parameter variation. X is the corres- 
ponding eigenvalue of the initial system. The derivatives of the eigenvalues 
to the system parameters are expressed in operational notation, r is the order 
of the derivative. 
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Truncating equation (8) after a desired derivative, the changed eigenvalue 
may be calculated approximately. Easy to handle is a linear formula with first 
derivatives. A more exact expression can be obtained if an additional quadratic 
term is taken into account. 


Linear Approximate Formula 

Taylor s expansion (eq. (8)), truncated after the first derivatives leads 
to the following linear approximate equation (9) 


8A 

A = A +3-5. 

1 


8A 


n no 9p 


Ap, + 


n 


1 3p, 


3A 


n 


K 


3A 

" ^no ^ 


n 


Apj^ (9) 


For determination of A respectively AA = A - A the first partial derivatives 
dA^/3p are required. n n no 

Using Oie standard subscript notation for partial derivatives 


3A 




n 


3M 


3C 


3k 


''n,k 3pj^ 3p^ -»k ■ ^'k ■ 3i^ -’k " 3i^ 

the differentiation of equation (3) with respect to p yields 

k 

In.k " ^ * *„£.k " £.k> 2„ ' ° 

( 11 ) 

We premultiply equation (1 1 ) by to obtain the scalar expression 


( 10 ) 


M . X^C . K) 5^^^ . il {2X„M . C) 3^ X___^ 


-*n ^-k * \ £-k * ^-k>3„ 


( 12 ) 


The first term in this equation is zero, it represents the left-hand eigenvalue 
problem, multiplied with q . Furthermore it is assumed that q and £ are 
normalized in a way to satiSty the relation 


£ {2A M + C) q = 1 

-n n- - -in 


(13) 


It follows that 

3A 


n,k 3p 


- £^ {A^ M, + A C,, + K„ } q 
n n — ’k n — ’k -’k _n 


(14) 
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Substitution of this result into equation (9) yields 


AX = (X -X ) = 
n n no 


- E + X C, + K, }q Ap, = E g Ap (15) 

-n nu’k no k - ’k -n *^k °nk *^k 


K 


with the influence coefficients 
Finally we introduce A-matrices 



K 

K 


AM = 

kJ, 

II 

II M 



K 

K 


AC = 

E 

= , E, 

C, 


k=l 

-k k=l 



K 

K 


AK = 

k=^l 

*!!k- A 


expressing 

the 

change of mass , 


k“^'k 


(16) 


the linear expression 


AX 


n LIN 


(X -X 


no^LIN 


{X^ AM + X AC + AK} q 
-n no — no— — ~n 


(17) 


This approximate equation is a good tool for calculation in many cases, pointing 
out the influence of parameter changes to the stability behavior of rotors. 

As above mentioned for application of the formula the modal parameters of the 
initial system are needed besides the parameter changes. Derivatives of the 
eigenvectors do not appear in equation (17). 


Influence Coefficients 


In equation (15) we have defined influence coefficients g^^ 


g. 


8X 

^ ni 

nk 8p, 


{X^ M,, + X C,, + K„ } q 
-n ncT— k nj^k -’k _n 


(18) 


This differential quotient, respectively differential sensitivity yields the 
influence on the eigenvalue X of an infinitesimal change of a particular 
parameter p^. It can be used as a meaningful approximation for finite parameter 
changes, when the considered modifications are relatively small. 


AXn 

^nk ^ Ap^ 


Aa Au) 

n ^ . n 


(19) 
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Influence coefficients 
of the initial system. 


are complex numbers. They depend on the modal parameters 


free 


( 


Changing the stiffness Ic of a spring with one end fixed and the other end 
see Table 1), we obtain with M,, = 0 and C, =0 

~ K — k 


i , = - a K,, q 

nk -n -*k 


( 20 ) 


The derivation leads to a unity main diagonal element K, , = 1 with all 

other elements cero and therefore the influence coefficient of the spring is 


®nk ^nk V “ (21) 


Table I contains further influence coefficients for the most important 
e ements in rotor dynamics like springs, dampers, masses, oil film bearings 

flnH CO o 1 c o 


Quadratic Approximate Formula 


into consideration the quadratic terms 
we obtain the following expression 


in Taylor’s expansion too, 


A 

n 


= A 

no 


K 

+ E 
k=l 





K 

E 

£=1 


2 

a A 


o 


( 22 ) 


in this improved formula the second derivatives A, are needed besides the 
first deriyatiyes and the given parameter changes. The determination of the 
second derivative is shown in the appendix. The expression is more complicated 
than the first derivative, but all quantities can be determined from the initial 


PARAMETER CHANGES AT SIMPLE ROTOR MODELS 


We applicate the method at simple rotor models to show the influence of 
particular parameter changes to the stability behavior. Naturally these simple 
examples can be solved also exactly, on the other hand they are suitable to 

emonstrate the approximate formulas. A comparison of the results with exact 
solutions IS possible. 
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Rigid Rotor with Flexible Supports and Clearance Excitation 


The first example is an elastically supported rigid rotor with mass m 
(fig. 3). The supports have the stiffnesses kj . in horizontal direction and 
in vertical direction. A clearance excitation force acts in the middle of the 
rotor. For instance such excitation forces appear in steam turbines. They result 
from the unsymmetrical fluid flow through the radial clearances at rotor and 
blading which appears according to the eccentricity between rotor and casing. 

The clearance excitation force acts rectangular to the displacement direction 
(fig. 3). 


r- -1 



The coefficient k depends on the power of a turbomachine, increasing with power. 

Investigating the translatory natural motion we employ the coordinates 
(j j and cj 2 • 


Equations of Motion, Eigenvalues and Eigenvectors 
The equations of motion for the rigid rotor with elastic supports are 


m 

o 


«s/ 

+ 

r 

L 

o 

m 


1 

to 


— 1 

CM 

CM 

M 

CM 

1 

1 


or with the definitions 



(24) 


2 

U) 

o 


3 


= 2kj j/m 

= k22/kj , 


= k/2k 


1 1 


T = cot () =d()/dT 

o 

stiffness ratio of support (25) 

dimensionless clearance excitation coefficient 


in dimensionless form 


1 0 


'll 



^ u 

0 1 


L'^2j 


~i 3 


1 

1 

1 

1 

no 

1 




(26) 
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The corresponding eigenvalue problem is (y = 


A/o) ) 
o 


1 + P 


with the characteristic equation 


+ (1+y) + y+B^ 


we obtain the solutions Z = \\ from 


= - ill + ■/ 0 -y)^-43^ 


respectively the four eigenvalues as a function of y and B. 

wi'rh 4 shows the real part and imaginary part of the essential eigenvalue 

n- 9 RW ^°"?tant. The system is unstable within the range of 

^ ~ . + 2p) . This range increases with increasing values of the 

clearance excitation coefficient B. Isotropy (y = 1) of the supports is the most 
disadvantageous case An increase of the anisotropy stabilizes the rotor system. 

For (1 23) > Y ^ (1 + 23) the real part of the eigenvalue is zero and the 

natural frequencies are split. 

amplitude ratio of the right-hand eigenvector corresponding to an 
tifopposite sign (27), the left-hand eigenvector hat 


_ i+y 


^2 ,.p2 

f, 6 


Parameter Variations 

To demonstrate the approximate formulas we choose an initial system with 
a^add^r investigate variations of y 3 aS 

for e?g«::?;e Zige^ a;e"“" -^-vectors the particular expreasiona 

“Variation of y with increment Ay 




o O fc 

o Ay q 


^ 2^2 
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-variation of 3 with increment A3 


Ay = - , 5^2^ 



(i^iq2 - 


-additional modal damping AD = c/2m03^ 



2AD o 

"'ll" 

Ay = - y [£j , 

0 2AD 



2y ^ 


The linear approximate formula points out the stabilizing effect of anisotropic 
supports (Ay positive in fig. 5b) . Figure 5a shows an increase of frequency 
caused by stiffening the system from y to y + Ay. 

For special cases the linear approximate formula has disadvantages, for 
instance in the isotropic case (y = 1 ) with horizontal tangent line. Then better 
results will be obtained with the quadratic formula (fig. 5b). 

The destabilizing effect of increasing clearance excitation coefficient 3 
and the stabilizing effect of additional damping AD are presented in figure 5c, 


Rigid Rotor in Journal Bearings 

Figure 6 shows a rotor of a test rig. Compared to the bearings the stiff- 
ness of the shaft is very high, therefore it can be regarded as a rigid rotor 
supported in two equal bearings (cylindrical bearings B/D = 0.8), rotating with 
the rotational speed Q. The rotor has a mass m = 72 kg, the distance between 
the bearings is 660 mm. The statical bearing load ^ ^ ^ equal to 353 N. Each 
journal has a diameter of 50 ram, the values for the radial clearance and for the 
oil viscosity are Ar = 210 ym, respectively 3.35 Ns/m^. 

The dynamic behavior of journal bearings can be characterized by four 
stiffness - and four damping coefficients and respectively by the non- 

dimensional quantities 


^ik 



Ar 

F 

stat 



'ik 


Ar Q 
F 

stat 


(31) 


They are functions of the Sommerfeld number. The pure translatory motion can be 
described by means of the coordinates q^ and q 2 * 


Equations of Motion, Eigenvalues and Eigenvectors 

The following equations of motion for a rigid rotor in two journal bearings 
describe the equilibrium of forces in the case of pure translatory motion and 
without external loads. 
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m 

o 



+2 






o 

m 


^2 



1 1 

^12 

21 

^22 



(32) 


Defining 


2 

( jO 

o 


g/Ar, T = CO t, w = fi/o) , ( ) = d ( ) / dT 

o o 


the equations of motion can be written in the nondimens ional form 


+ 1 
w 


1 1 


21 


12 


22 


The statement for the natural motions 




1 1 


'21 


12 


'22 


= 0 
(33) 


Xt 

q. = q. e 
1 1 


yx 

e , y 


= X/co 


(34) 


yields the eigenvalue problem 

^11 2 
Yi 1 + ^7- u + y 


3 


12 


1 1 w 


3. 


^ 21 
^21 ^ 


^12 " w 

3. 

'22 w 


^ 22 2 

Yoo + TT- y + y 


= 0 


(35) 


The four eigenvalues y^ are i. g. conjugate-complex quantities. If the eigen- 
values are known, the corresponding eigenvectors can be calculated from equation 
(35). 


Figure 7 shows the eigenvalue that causes the instability in dependence 
the dimensionless rotational speed w. Also the eigenvalue that belongs to a 
rotatory motion is plotted, but it has no meaning for the stability behavior of 
the rotor. The real part of the eigenvalue will come to zero at w = 2.66. 

At this rotational speed the stability threshold is reached. The Sommerfeld 
number here has the value So = 0.563. The corresponding circular natural fre- 
quency ^ 1 ^^ is 1.435. 
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Influence of Parameter Changes at the Stability Threshold 

We start from the stability threshold and we investigate how particular 
parameters effect the stability behavior. Therefore we take the linear approxi- 
mate formula (17). First we ask about the influence of the single coefficients 
Yik and 3»j^/w. It can be judged by the influence coefficients in Table 1. Figure 
8 shows tfte influence coefficients in the complex plane. For example, one can see 
that an increase of y ^ is labilizing or that an increase of the damping coeffi- 
cients stabilizing the rotor motion. 

If we change the rotational speed, the Soimnerfeld number and the stiffness 
and damping coefficients also will change. Figure 9 shows the variation of the 
eigenvalue in dependence of the rotational speed calculated by the approximate 
formula. The formula figures the tangent line to the eigenvalue curve at the 
expansion point. An increase of the rotational speed causes instability. 

Finally in figure 9 the influence of the oil viscosity is shown. If we 
increase the oil temperature from 38°C to 44®C, the oil viscosity will decrease. 
The real part of the eigenvalues becomes negative and there is a stabilizing 
effect . 


ELASTIC ROTOR WITH FLEXIBLE SUPPORTS AND SEAL FORCES 


A more complicated rotor-system, an elastic turbopump rotor with flexible 
supports, an impeller and two plain seals is shown in figure 10. Between the 
bearings where the basic shaft diameter is 64 ram, the shaft carries the impeller 
mass of 55 kg. Two plain seals are mounted besides the impeller with a seal 
length of 40 mm and a radial clearance of 0.3 mm. The dynamic characteristics of 
the seals are given by the stiffness, damping and inertia coefficients K, k, C, 
c, m* (fig. 10). The shaft is supported in two identical bearings, each one having 
a horizontal stiffness k = 1.0 107 N/m and a vertical stiffness V ,22 
The distance between the oearings is equal to 1.2 m. 

First of all eigenvalues and eigenvectors were calculated for the described 
original system in a speed range from 2000 rpm to 8000 rpm. In figure 1 1 two 
damping constants (decay constants) are plotted versus rotational speed. One of 
the damping constants crosses the zero axis, pointing out the instability onset 
speed of 5200 rpm, which is above the operating speed. 

In order to improve the stability behavior, some parameter changes were 
investigated, starting from the initial system with a rotational speed of 5000 rpm 
Figure 11 shows the results, especially the influence of seal stiffness and dam- 
ping coefficients. All coefficients were changed 20 per cent. An increase of K 
and C stabilizes; an increase of k has a destabilizing effect. 

Furthermore it is shown that a positive change of k ^2 direction to 
isotropic bearings is disadvantageous. 
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VARIATION OF THE MASS OF A TEST RIG ROTOR 
WITH UNSYMMETRIC FLEXIBLE SHAFT 


The test rig rotor shown in figure 12 consists of a flexible shaft with 
rectangular cross section 8 x 12 iran . It is running in two ball bearings which 
are comparatively stiff. In the middle between the bearings a disk is supported 
having a mass of 0.89 kg. The nonrotating shaft has two different natural fre- 
quencies 24.6 Hz and 35.1 Hz concerning to the different stiffnesses in two 
rectangular planes. 

It is well known that in rotors with unequal moments of inertia different 
dynamic effects may occur, for instance vibrations caused by the rotor weight 
and unstable vibrations (ref. 6) . Concerning the stability behavior in this case 
a variation of the disk mass was investigated by calculations (exact results and 
approximate formulas) and for control by measurements at the test rig. Theore- 
tical results for the eigenvalues of the undamped shaft are plotted in figures 
13a and 13b versus the rotational speed and the mass of the disk. The eigenvalues 
defined for the rotating coordinate system are either pure imaginary (fig 13a) 
or pure positive real (fig. 13b). The last one occur in a rotational speed range 
between the two natural frequencies of the nonrotating shaft. That means that 
the rotor motion is unstable in this range. 

In order to show the utility of the approximate formula the change of 
the eigenvalues (rotating system) caused by a positive 25 per cent variation of 
the disk mass was investigated for the three different running speeds 1200 rpm, 
1800 rpm and 2400 rpm. The three figures 14a, 14b, I4c show the variation of 
the natural frequency, respectively the damping constant, in dependence of the 
mass. The mass of the original system is 0.89 kg, the changed mass is 1 . 1 kg. 

The three diagrams contain the exact, the linear and the quadratic approximate 
solution. Furthermore the results from measurements are plotted in figures 14a 
and 14c. There is a good agreement of the results. For 1200 rpm the mass change 
decreases the frequency (destabilizing effect), for 2400 rpm the mass change 
increases the frequency (stabilizing effect). A stabilizing effect is also 
given at the running speed 1800 rpm, where the decay constant a is lowered by 
the mass change. 
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APPENDIX 


Calculation of the Second Derivatives of the Eigenvalues 

First we have to express the right- and lefthand eigenvalue problem 
3) » (eq. 7) in the following form: 

(A 1) 
(A 2) 


(A - X B) X = 0 
— n- — n 

W - - 0 


where the matrices A, B and the vectors x and y are defined as 

— n ^T\ 

M 0 


A = 


0 -K 


“n '^n 
0 M 


B = 


M C 


X = 





X Z 

= 

n~n 

’ ^n = 

n-n 


J 

£ 

n 


The vectors x^^ and satisfy the biorthogonality relation 


(A 3) 


(A 4) 


y B X =6 
— n — — m nm 


(A 5) 


where 6^ is the Kronecker delta. The first derivative of X with respect to the 
parameter p^ has the form 


X . 
n.J 



B X ) 
n 



(A, - X B,.) X 

- J n— j -n 


(A 6) 


For the calculation of the second derivatives of X , we need the derivatives of 

the eigenvectors _x and _y . The vectors x . and*^y . can be represented as a 

linear combination of x and v : 

— n — n 


X 

~n.J 



c . . X. 
njL -1 


2N 

.1 d . . y. 
1=1 nji -^1 


(A 7) 


Using the equations (Al), 
the following expressions 


A2), (A5) and (A7), Plant 

for the coefficients c . . 

nji 


and Huseyin (ref. 3) obtain 

and d . . 
nji 
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C. +d. ~ ~ 2 . 2£n 

njn njn -^n ~ j 


c . . 
nji 


y. (A, .-A B, . )x 
^ 1 — Vi n-*j -n 


X - X. 

n 1 


d .. 
nji 


y^(A, .-X B,.)x. 
-n - J n-*j -1 


X - X, 

n 1 


I for tii4 i 


(A 8) 


In order to get the second derivative of X with respect to p. and p^^ we must 

differentiate eq. (A6) with respect to p^^.^With eq. (Al), (A2) , (A5) and 

(A7) we obtain for X . p (ref. 3) 

' n . 1 it 




n 




X 

— n 


+ (c . + d . ) X n + (c p + d„p„) X . 

njn njn^ niln nx.n n,j 


2N 


+ .Z, (X - X.) [d .. c p. + d p. c ..] 

1=1 ^ n X nji nici nici nji 

i+n 


(A 9) 


where c and d .. can be calculated by eq. (A8) . 
nj 1 nj 1 

If we go back to our initial problem with the matrices M, C and K and 
the eigenvectors ^ and we get by substituting A, x and ^ for the second 
derivative 


3X 


^ . _J1_ — = - (X^ M, .p + X C,.p K.^o) q„ 

n,jii 3p*8pn "n n — jil n jA j>c n 

J ^ 

- ( * £.j) a„ > 

- i * C,j) ^ 


ii, 

i+n 




where G . means 
-nj 


(A 10) 
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G . 

-nj 


r M, . + A C, . + K, . 
n - J n-’j -’j 


(All) 


The second derivative with respect to a single parameter p. can be obtained 
by replacing the subscript £ in eq. (AlO) by the subscript*^ j • 
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Figure 6. - Rigid rotor in journal 
bearings. 


Figure 7. - Eigenvalues of rigid 
rotor in journal bearings. 
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Figure 8. - Influence of bearing 
coefficients to stability 
behav ior . 


Figure 9. - Variation of rotational 
speed and oil viscosity. 
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Figure 10. - Elastic turbopump rotor. 



Figure 11. - Influence of seal stiffness and damping coefficients 
to stability behavior of turbopump rotor. 
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